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NOMENCLATURE 


English Symbols 


A - Area (eee) 

Agr - Drum Downcomer Total Cross Sectional Area (£t7) 

Aas - Header Downcomer Total Cross Sectional Area (£t*) 

AMB - Main Bank Heat Transfer Surface Area (ft?) 

A. - Main Bank Riser Total Cross Sectional Area GRE 

A. - Screen Riser Total Cross Sectional Area ge) 

Be - Becker's Factor 

SC - A General Constant 

Cy - Specific Heat of Liquid (Btu/lbm/°F) 

Cn - Specific Heat of Tube Metal (Btu/lbm/°F) 

cS - Specific Heat (Btu/lbm/°F) 

Cel - Specific Heat of Flue Gas in the Furnace 
(Btu/lbm/ °F) 

Cela - Specific Heat of Flue Gas at the Furnace 
Exit (Btu/lbm/ °F) 

CR> - Specific Heat of Flue Gas at Superheater Exit 
(Btu/lbm/ °F) 

CR3 - Specific Heat of Flue Gas at Main Bank Exit 
(Btu/lbm/ °F) 

as - Constant used in Rohsenow Correlation 

C.-C) 4 - Constants Used in Various Equations 

d - Diameter (ft) 

Dar - Sum of the Diameters of the Drum Downcomers 
(EE) 

Das = a the Diameters of the Screen Downcomers 





ED1 


ED2 


ER1 


ER2 


ER3 


ER4 


wr 


EW1 


EW2 


EW3 


EW4 


Sum of the Diameters of the Main Bank Tubes 
CE) 


Sum of the Diameters of the Screen Tubes (ft) 
Total Heat Input to the Furnace (Btu/sec) 


Heat Transferred from the Desuperheated 
Steam to the Tubes (Btu/sec) 


Heat Transferred from the Superheater Tubes 
to the Steam (Btu/sec) 


Heat Transfer Rate from the Flue Gas to the 
Screen Tubes (Btu/sec) 


Heat Transfer Rate from the Flue Gas to the 
Superheater (Btu/sec) 


Heat Transfer Rate from the Flue Gas to the 
Main Bank (Btu/sec) 


Heat Transfer Rate from the Flue Gas to the 
Economizer (Btu/sec) 


Heat Transfer Rate to the Riser Liquid 
(Btu/sec) 


Heat Transfer Rate from the Economizer Tubes 
to the Feed Water (Btu/sec) 


Heat Transfer Rate from the Desuperheater 
Tubes to the Drum (Btu/sec) 


Heat Transfer Rate to the Screen Riser Liquid 
(Btu/sec) 


Heat Transfer Rate to the Main Bank Riser 
Liquid (Btu/sec) 


Friction Factor 

Drum Downcomer Friction Factor 

Heater Downcomer Friction Factor 

Main Bank Riser Friction Factor 

Screen Riser Friction Factor 
Acceleration of Gravity - 32.2 ft /sec* 


Mass Flow Rate (lbm/sec) 
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GB 
GDl 
GD2 


GD3 


GD4 
GL 
GR 
GS 
Gr 
Gw 
GP 


Gwr 


Gws 


HD1 
Hig 
Ho 
HNB 
HNBR 


HNBS 


Fuel Flow Rate (lbm/sec) 
Total Steam Flow Rate from Boiler (lbm/sec) 
Desuperheated Steam Flow Rate (lbm/sec) 


Steam Condensing from Vapor to Liguid in the 
Steam Drum (lbm/sec) 


Superheated Steam Flow (lbm/sec) 

Air Mass Flow Rate (lbm/sec) 

Total Flue Gas Flow Rate (lbm/sec) 

Screen Riser Circulation Flow Rate (lbm/sec) 
Riser Circulation Flow Rate (lbm/sec) 
Downcomer Circulation Flow Rate (lbm/sec) 
Main Bank Riser Mass Flow Rate (lbm/sec) 


Drum Downcomer Circulation Flow Rate 
(lbm/ sec) 


Header Downcomer Circulation Flow Rate 
(lbm/sec) 


Enthalpy of Saturated Steam (Btu/lbm) 
Heat of Vaporization (Btu/lbm) 

Heating Value of Fuel Oil (Btu/lbm) 
Non-boiling Height (ft) 

Main Bank Riser Non-boiling Height (ft) 
Screen Riser Non-boiling height (ft) 
Enthalpy of Saturated Water (Btu/lbm) 
Enthalpy of Downcomer Liquid (Btu/lbm) 


Enthalpy of Feed at Economizer Inlet 
(Btu/1lbm) 


Enthalpy of Feed at Economizer Outlet 
(Btu/lbm) 


Enthalpy of Water in Lower Drum (Btu/lbm) 


Tey 





Hxr 


Hxs 


k 

K 

Ks 
KR1-KR4 
R 

L 

LB 
Ldr 
Lds 
Ldlr 
Ldls 
tr 


Lrl 


Lsl 


LMTD 


Effective Height of Main Bank/Drum 
Downcomer Circuit (ft) 


Effective Height of Screen/Header Downcomer 
Civeume (ff) 


Thermal Conductivity (Btu/hr-ft-°F) 

a general constant 

Equivalent Sand Roughness (in) 

Constants Used in the Gas Heat Transfer Equations 
Length (ft) 

Tube Length (ft) 

Boiling Length (ft) 

Total Length of Drum Downcomer Tubes (ft) 
Total Length of Header Downcomer Tubes (ft) 
Avg. Length of Drum Downcomer Tubes (ft) 
Avg. Length of Header Downcomer Tubes (ft) 
Total Length of Main Bank Tubes (ft) 

Avg. Length of Main Bank Tubes (ft) 

Total Length of Screen Tubes (ft) 

Avg. Length of Screen Tubes (ft) 

Log Mean Temperature Difference (°F) 

Mass (lbm) 

Mass of Screen Tube Metal (l1bm) 

Mass of Superheater Tube Metal (lbm) 

Mass of Riser Tube Metal (1bm) 

Mass of Economizer Tube Metal (1bm) 

Mass of Desuperheater Tube Metal (1bm) 
Mass of Flue Gas in Furnace (lbm) 


Mass of Water in Steam Drum (l1bm) 


We 





PD 


PD1L 


PD2 
PD2 
Pr 
Psat 
PW1 


PW2 


Ral 


Ra2 


Re 
RW 


RS 


RO 


Tsat 
TD 


TDL 


TD2 


Pressure (psia) 
Pressure of Steam Volume in Steam Drum (psi) 


Saturation Pressure Corresponding to Enthalpy 
of Water in the Steam Drum (psi) 


Pressure at Superheater Outlet (psia) 
Pressure at DeSuperheater Outlet (psia) 
Prandtl Number 

Saturation Pressure (psia) 

Pressure in Water Drum (psi) 

Pressure in Screen Header (psi) 

Heat Transfer Rate (Btu/sec) 

Tube Radius (in) 


Average Density in the Main Bank Risers 
(Lbm/£t3) 


Average Density in the Screen Risers 
(lbm/f£t3) 


Density of Saturated Steam (lbm/£t?) 
Reynolds Number 
Density of Saturated Water tiene) 


Density of Two-Phase Mixture at Screen 
Outlet (lbm/ft?) 


Density of Two3Phase Mixture at Main Bank 
Outlet (lbm/ft°” ) 


Time (sec) 

Temperature (°F) 

Saturation Temperature (°F) 
Temperature of Steam Drum Vapor (°F) 


Temperature of Saturated Steam Corresponding 
to Drum Enthalpy (°F) 


Superheater Outlet Temperature (°F) 


LS 





TD3 
TM1L 
TM2 
TM3 
TM4 
TM5 
TRL 
TR2 
TR3 
TR4 
TR5 


TRa2 


TRa3 


TRa4 


VB 
VDL 
VENT 
VNDR 
VR 


VS 


Desuperheater Outlet Temperature (°F) 

Screen Metal Temperature (°F) 

Superheater Metal Temperature (°F) 

Main Bank Metal Temperature (°F) 

Economizer Metal Temperature (°F) 

Desuperheater Metal Temperature (°F) 
Temperature of Flue Gas in Furnace (°F) 
Temperature of Flue Gas Leaving Furnace (°F) 
Temperature of Flue Gas Leaving Superheater (°F) 
Temperature of Flue Gas Leaving Main Bank (°F) 
Temperature of Flue Gas Leaving Economizer (°F) 


Average Flue Gas Temperature Across the 
Superheater (°F) 


Average Flue Gas Temperature Across the 
Main Bank (°F) 


Average Flue Gas Temperature Across the 
Economizer (°F) 


Feed Inlet Temperature (°F) 

Feed Outlet Temperature (°F) 
Ambient Temperature 

Bubble Rise Velocity (ft/sec) 

Mean Velocity (ft/sec) 

Volume feo) 

Boiling Volume (£7) 

Volume of Steam in Steam Drum (£>) 
Percent Throttle Valve Opening 


Water Drum Volume (ft?) 


Main Riser Bank Fluid Volume (€2>) 


Screen Riser Bank Fluid Volume (£t?) 
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cy = 
X = 
XD = 
XS = 
Xe = 
Z1 = 
Z2 = 
Z3 - 
ZA - 
Z5 - 
Z6 = 
Z7 = 
Z8 - 


Specific Volume (€t?/1bm) 

Quality (percent) 

Quality at Main Bank Riser Outlet (percent) 
Quality at Screen Riser Outlet (percent) 
Quality at Riser Exit (percent) 

Drum Downcomer Bend Loss Factor 

Drum Downcomer Exit Loss Factor 

Main Bank Riser Bend Loss Factor 

Main Bank Riser Exit Loss Factor 

Screen Riser Bend Loss Factor 

Screen Riser Exit Loss Factor 

Header Downcomer Bend Loss Factor 


Header Downcomer Exit Loss Factor 


Greek Letters 


u = 
0 = 
oO = 
°pal 7 
°Da2 

Subscripts 
a _ 
ie, aes 
M,m - 


Viscosity (lbm/ft-sec) 

Density Gis 

Surface Tension (1bf/ft) 

Average Density in Superheater (1bm/£t7) 


Average Density in Desuperheater (lbm/£t?) 


Average 
Liquid 
Vapor 


Metal 


Mts) 





ar 


ds 


Far Field Average 
Main Bank 

Screen 

Drum Downcomer 
Header Downcomer 
Steam 

Flue Gas 


Water 


IE 





I. INTRODUCTION 


The increasing attractiveness and efficiency of computer- 
aided design techniques have established a need for 
reasonably accurate computer models of all types of 
engineering systems. The design of multivariable control- 
lers, in particular, requires the availability of state 
information about a system which is not available from 
conventional transfer function simulations. A number of 
computer models of marine boilers have been proposed 
{1,2,3,4]. All of these efforts share one or more of the 


following features: 


a. Lumped parameters are used. 


b. Bubble formation effects with respect to shrink and 
swell are not modeled. 


c. The furnace screen and main riser bank are lumped 
together as are their respective downcomers [1,3,4]. 


d. The model is linearized using a Taylor expansion 
Or perturbation technique [3]. 


Use of lumped parameters is dictated by the need to maintain 
computational efficiency, ero anllenoey when the model is 
used in control system design and optimization schemes. 

In the past, effects due to vapor formation in the drum 

have not been included due to the lack of a suitable model. 
A simple model based on the Rohsenow correlation [5] was 
developed for use in the current simulation. However, its 


contribution to shrink and swell effects was negligible and 
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it was eliminated from the equation set. Regarding the 
tube banks, Haraguchi and AkaSara [2] split out the various 
riser and downcomer loops, including intermediate headers. 
Considering the general level of uncertainty inherent in 
the correlations used in this model, this appears to be an 
excessive degree of complication. A compromise position 
was taken in the present effort in that the screen and main 
riser bank and their associated downcomers were modeled 
separately, but intermediate headers within the screen 
circuit were not addressed individually since their configura- 
tion varies widely from one boiler to another. 

The objective of this thesis was to develop a general, 
non-linear model of a Navy boiler for use in simulations 
covering the full operating range of the propulsion plant. 
However, Since a linearized version 1s normally needed for 
multivariable controller design, CENTRL, a computer code 
developed locally for use with CSMP models, could be used 
to generate linearized state space representations described 


by the matrix equations: 


Ix} = [a] fax} + [8] fouy 
Fayt= Le] fax} + [vo] pou} 


This approach would produce linearized representations 


covering the entire boiler operating range in a piecewise 
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fashion, rather than attempting to extend the range of 
applicability of a single state space or linear model. 
Since this author found the notation used in Reference 


[1] particularly readable, it was adopted for use in this 


document. 


aL, 





II. DESIGN CONSIDERATIONS 


A. D-TYPE BOILER 

The single furnace (D-Type), uncontrolled-superheat 
boiler was chosen as the basis for the model since it is 
the predominant design in the current U.S. Navy inventory. 
The nominal saeratiine pressures range from 600 to 1200 psi, 
and all designs are conceptually similar. A basic schematic 
is shown in Figure 1. These boilers are typified by 
relatively fast response times for use in high speed 
maneuvering vessels. This puts heavy emphasis on properly 
designed control systems. Manual control of these boilers 
is essentially infeasible except under the most favorable 
circumstances, e.g. steady steaming at low speeds. Current 
conventional control systems used by the Navy are of the 
Same generic type with relatively minor customization except 
in hardware from one installation to the next. These 
designs have been based on analog models intended to provide 
reasonable dynamic responses for pressure, water level, and 
steam flow which are the only state parameters measured by 
conventional controllers. Although these control systems 
have in general proven satisfactory, there is some indica- 
tion [6] that multivariable control systems can effect 
Significant improvement in both response time and demands 


on auxiliary equipment. 
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A typical boiler of this type, the Combustion Engineering 
Type V-2M is shown in Figure 2. The main components of the 
boiler are the furnace which is surrounded on all sides 
except the floor by screen tubes, the superheater, main 
riser banks, economizer, water drum, steam drum, desuper- 
heater, downcomers, and screen headers. Navy boilers 
do not include air heaters and the only air preheating 
which occurs is accomplished by passing the incoming air 
over the exhaust stack liner. This procedure generally 
produces an air temperature of about 100-130°F for an 
ambient of 70°F. In contrast to the boilers used as the 
basis for previous models, the desuperheater is located 
in the water drum rather than in the steam drum. All of 
the steam generated passes through the superheater prior 
to being split out as main superheated steam and auxiliary 
desuperheated steam. 

All of these boilers are of oe Nacura se Eculation 
type. The “throttle valve" was considered to be in the 


Superheated main steam line. 


B. MODEL FLEXIBILITY 

The model was designed to be as flexible as possible 
to permit its relatively simple adaptation to a particular 
boller. Since the boiler technical manual or equivalent, 
the steam and gas tables, and standard engineering hand- 
books may be presumed to be available, the model was designed 


to utilize data only from these sources for implementation. 
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To convert this input data into the form required by the 
model, an interface program was developed to adjust’ 
parameters and determine initial conditions at the operating 
point for which the input data were extracted. This 
approach was taken to avoid the necessity of redesigning 

the model for each operating range or specific boiler 


design being considered. 


C. ASSUMPTIONS 
To permit tractable computation, a number of simplifying 
assumptions were made in the formulation of the model: 

1. Feed water entering the drum mixes thoroughly with 
the drum liquid. 

2. No steam generation occurs in the steam drum except 
through mass transfer at the water-vapor interface. 

3. There is a uniform distribution of heat flux in the 
generating tubes. 

4. The mixture of steam and water in the generating 
tubes is homogeneous. 

5. Quality varies linearly in the generating tubes. 

6. Heat losses to the ambient except through the 
exhaust gases are neglected. 

7. All downcomers and risers in a particular circuit 
have the same length as the average length, and the 
shapes and diameters are equal. Also the number of 
tubes is the same as in the actual boiler. 


8. The slip ratio is 1.0 (no slip). 
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' 


Since the feed enters through a long, horizontal feed 
pipe, and the drum liquid iS in a state of turbulence, 
the first assumption is well justified. Although the 
second assumption is not strictly correct during transients 
when there is a Significant change in drum pressure, and 
hence saturation temperature, the effect of steam generated 
at the drum wall-liquid interface appears to be negligible 
based on an analysis of this factor in connection with 
this study. The assumption of uniform heat flux is 
recommended by Reference [7] aS a reasonable alternative 
to a stepwise integration over the tube height. The 
assumption of a linear change in quality is a consequence 
of the assumption of a linear heat load. For a well 
designed boiler, heat losses through the casing are small 
compared to the total energy generated. The effects of 
any such losses will be reflected in small errors in the 
computed heat capacities of the gas and in the heat transfer 
coefficients. The assumptions regarding tube geometry are 
made primarily for computational simplification. The 
last assumption regarding slip is not an essential one, 
and higher slip ratios could easily be considered if 
necessary to the analysis of a particular boiler design. 
In addition to these explicitly stated assumptions, those 
inherent in the correlations employed in the model formula- 


tion are also considered to have been made. 


Zs 





IIIT. MODEL ANALYSIS 


A. GENERAL 
The model was envisioned as an assemblage of control 
volumes for which the rates of energy transfer are 
described by appropriate equations to obtain the rate of 
change of energy storage. The resulting differential 
equations were integrated to obtain the "state" of the 
system. In this way, the states can be chosen to have 
specific physical interpretations. The system is shown 
diagrammatically in Figure 3. 
The following control volumes were considered in the 

model formulation: 

1. Riser tube walls 

2. Superheater tube walls 

3. Economizer mass 

4. Desuperheater tube walls 

5. Riser fluid kinetic energy) 

6. Riser fluid (thermal energy) 

7. Steam drum water volume 

8. Steam drum steam volume 

9. Water drum 


10. Furnace 


B. FLUE GAS AND COMBUSTION SYSTEM 
The flue gas and combustion system is shown schematically 


in Figure 4. Heat is carried into the furnace through the 
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sensible heat of the air and fuel, and the heat released 
in combustion. Since the fuel-air ratio is relatively 
constant, and the fuel and air have relatively constant 
inlet temperatures, the sensible heats were lumped with 
the heating value of the fuel. The mass flow rate into 


the furnace is given by 

Gea=Gg + Gy 
and the total heat input is 

Ea = Gaxo 


The heat transfer to the furnace screen is primarily through 


radiation, 1.e., 


+ 460\" Ti 5 ( 
Ea, 3 ae ae! = ae) 


/00 1/00 
The heat balance for the mass of flue gas in the furnace is 


< (Med Cea Tes ) = Eg -Gai - Ga Ce, Tey 


where the last term represents the heat carried off by the 
Flue gas. Heat transfer to the remaining tube banks is 
primarily through convection. The heat given up by the 


flue gas 1S given by the general relation 


$335 = M Ce La: 
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and the heat transfer to the tubes may be estimated with 


the Grimson Correlation [8] 


Ptobe = 76 Re F Bp (7% - Tm ) 


Since the physical properties in the above equation do 
not change significantly over the expected range of 
temperature variations, the correlation is used in the 


form 
Fire = KGa CT, - Tm) 


Since the tube geometry 1S not known in advance, ann = 0.6 
was used as a reasonable average for most geometries found 


in practice. The equations for the tube banks were: 


Ll. Superheater 


Ge Cerio (Taz - 181) 
'<"F5 Ga (Trar —Tmz) 


VQaz2 (TR2 + Teas) [2.0 


ERx 
ER2r 
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2. Main Bank 


Eas = Gpelerr (Taz —7R4) 
Fes = Ke, Ge (Gress — Tm3) 
VRAZ= (Taz + Tay) /2-0 


3. Economizer 


Ery = Ge Cez (Tey -7es) 
Eeq4 = KRy Ge" (Trad ~Tiny) 
‘TRay = (Try + Tas) [ a.0 


The arithmetic average temperature for the gas was used in 
the driving potential AT rather than the log mean 
temperature differences (LMTD), because a comparison of 
the results obtained for the two methods indicated the 
differences were slight for the temperatures involved and 


the simpler average was more computationally efficient. 


C. WATER-STEAM CIRCULATION SYSTEM 
As is shown in Figure 5, the basic circulation system 
consists of the steam drum, water drum or header, and the 


riser and downcomer tubes. The driving force for circulation 


a 





is the difference in density between the riser and down- 
comer loops. The pressure loss around the circuit balances 


this driving force. For the downcomers, the force balance is 





2 
+z, +Z, +1) Gwr 


Pwa - Poa = Rw 9 Axe - (Fdy - A Ady? Rw 


for the drum downcomers, and 


G 2 
Pw 2 - Psi 2 Rw 4 Hxs oa ( FJ. Wg +1,+Z bo ra ie 


for the Screen downcomer. 


For the risers, the equations are: 


_>? = (¢ be Go Rad Gur” * 
Pua ~ Poa = (Fes + 1); 42 po? tye © Li ogaeer ae 





+ Rarg Hart Ras So 2 Cas. 


2 2 
Pusan — Poa =( Fe 2 Ree je hax +k, _ Gus +h Gs 


Re iy Py ad) al 2 
lsa of (Gs 
+ far 9 Has + laa 2S je (F) 


Rs 


for the screen and main bank, respectively. The friction 


factors are calculated from [9] 


5 [ete 
: ] C174 — 2 leg Fx) 
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The friction factors for the risers are multiplied by 


Becker's factor [10] to account for two-phase effects 


96 


Be=L + ieo78 28 (5) 
Pos 


The relationship between riser and downcomer mass flow 


rates is given by the continuity equation 
oS ier gee 
de 7°" V 


The average density in the risers is found from 


\ 
Wee = | { PUD AL rR, Hoa 
| “Hwa 


From the assumption that the quality varies linearly, the 


density may be obtained as a function of 2 from 


- L 
YD= Remedi HME AS 


Regrouping the terms we get 


uh 


QQ) = ———_—_.- a 3 
P (Ve _ ac we V5q ) + He Vis d 


Zo 





This expression can be integrated in closed form to yield 


Xx 
Wak cde, | et = 
Ke Vig EY + Vt — 4 Ane Veg 


Os) 





Combining terms and simplifying finally yields 


L 
Le e Vtg 
| pI = Za I sey 
H 


NB 
From this result, ee ls obtained as 


ay = aus _ +8 ln eae My 
% [L z eq + fy Ng 


An energy balance on the riser tube banks yields the 


quality 





S ¢ 
S| Sav ( Hus 7 ee tg) Ve 


= (Ear Awe a Ce (wd + Xe Hey ) + Ewr - Cw (Hws ~Hwr) 
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where VB is the boiling volume, and the last term on the 
right accounts for heat added in the non-boiling height 
to raise the enthalpy to saturated conditions. 


The non-boiling height is obtained from 





H ne =L om 


where qe is the sensible heat added to the liquid in the 
non-boiling height, and Gy 1s the total heat input. The 
ratio Go/ 4, 1S given by 


Rwi 7 Hwr 
( Hut + Xe Htg ) - Huv 


The boiling volume is related to the non-boiling height 


by 


Ve = L— Hae V 


With substitution of the appropriate parameters, the above 
equations were applied to both circulation loops. 
The mass balance of the water in the steam drum is 


d Mwt 


aE = Go (i= %e) * Gs (1-Xs) + Gos + Gwe —Gwr ~Gus 


ak 





and the energy balance is given by 


7 (Hw1 Mws | = Go (\-x.) Hwi + Gs (i—Xs) Hwa 
+ Gp3 H tq + Gw2 Hwy — Gur Hw — Gws |tw2 


Since all the incoming feed must be passed through the 
riser tubes to be heated, an energy balance on the 


downcomers gives 


Gwri Hwy + C Gur + Gws - Gw2 ) Rw 
a C Gur + Gws ) Hw2 


Gwr (Hwy —Hwi) + CGuwr + Gus) Hwa 
= C Gur i Gws ) Hw2r 


Solving for HW2, the downcomer enthalpy, we obtain 


Gw2 
tw. = -— 


= Culp b Gas ( Hw — Hwa ) = Hw 
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D. WATER-STEAM HEAT TRANSFER 
l. Risers 

The heat transferred to the riser tubers is used to 
bring the incoming flow to saturation conditions in the 
non-boiling height, and then generate steam in the remain- 
ing portion of the tube. To permit a tractable equation 
set, however, the tube walls were assumed to be at a uniform 
temperature, and all the heat was assumed to have been 
transferred under fully developed flow boiling conditions. 
Since the non=boliling height for the main bank was found to 
be small, the assumption was reasonable for that circuit. 
The non-boiling heights were higher in the screen circuit; 
however, the error in temperature introduced was less than 
three sapere which was deemed an acceptable level. For 


fully developed flow boiling, Tong [11] recommends 


4/3 
— Psat 3 
F ~ 1 1BR x 10% A Vset 


Therefore the heat transfer to the risers was given by 


4/3 3 
Ew3 = Kw3 Pow Cit -To1) 
Ew = Kwa Por 2 CT m3 Pa ye 


The tube metal temperatures were obtained from an 


energy balance, 


32 





Ma. on oie Sregi + Gw3 


for the screen, and for the main bank 


Van STs = Eg: — Ewa 





Cae Desuperheater 


The Dittus-Boelter correlation [12] was used to 


determine the heat transfer from the steam to the tubes, 1l.e., 


*& 
Eva > Kor Gor Cumrp) 


The log mean temperature difference was 


(Tox = To3 ) 


in oe i og 


(Tb3 — Tos) 


LMTD = 


The energy given up by the steam was obtained from 


Gor = Cor Gor CT. -T3 ) 


The heat transfer to the water in the drum was 


Ewt = Kwrn CIms - Tw ) 


34 





The metal temperature was obtained from the energy balance 


Eo. — Ewr 





il 


Mms Cm ws 


3. Economizer 
The Dittus~Boelter correlation gives the heat 


transfer from the tubes to the water as 


8 
Ewi = Kw GS w2 LMTSD 


Twa m w3 


LMTDH = 
In {a -Tw 3 


Tey —Tw4 
The heat absorbed by the water was 
Ewa = Gua Cw (Ty — Tw; ) 


4. SsUPperheater 


The heat transfer to the steam was given by 


.8 
Eht = Kyt Gow CTy2 = Tha) 


and the heat transferred from the metal in the tubes 


Was 
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*g 
Eonar > Kos Goy LMT 


Tou = Tot 


Tm -Tox~ 


LMrTbD = 


The metal temperature was found from 





Mma Cm Tm = Eerw- Eo 


E. PRESSURE DROP CALCULATIONS 


For turbulent tube flow, the pressure drop is given 





by 
L Det 
Ar = \ /» Gay ® Ye 
Now 
CS 
Um = A Oay ; 
SO 
2 
AP = ny eee ae 


A? Pav 


Lumping the constants we get 


4 
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z 
AP=zC S ; 


Therefore the pressure drops in the superheater and 


desuperheater were given by 


z 
Pot - Pow = C4 a 
ad 


and 


2 
“Por - Pr a Ca Goa 
Far 


The average densities were found from the ideal gas law 
under the assumption that the temperature variation can. 

be neglected, 1.e., density is a function of pressure only. 
This assumption is motivated solely by the need for compu- 
tational efficiency. It is reasonable for the desuperheater 
Since the temperature excursions there are small. Although 
errors are larger for the superheater density values, they 
are still within accuracy overall limits for the model. 


The flow through the throttle value was given by 


Goy = Ca Pox VENT 
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F. EQUATIONS OF STATE 

The following equations, obtained from Reference [13], 
were used to determine state properties. Within the range 
of pressures for which the model is valid (300-1500 psi), 
the results were within a few percent of tabular values in 


the worst case. 


“Psat 


Tsar 


Exp] C In Hsat — 4-4u1068) /.2bus% | 


ExP{ rzisa In Psat + 4:77:23 f 


‘mae 
Ts = 63-8 - O.OIMBL Isot + INZRE-S 'sat 


3 
— 6@18GE-8 Teat 
Veg = SA4/ Peat — O4 
2 
hg = 422.13 — 405i Peat + Utne -4 Peat 


3 
— 4.219 E-8 “Psat 


38 





IV. DETERMINATION OF CONSTANTS AND INITIAL CONDITIONS 


A. GENERAL 

For the model to be useful, it should be possible to 
obtain the necessary constants and initial conditions from 
the data supplied in the boiler Technical Manual augmented 
with reasonable engineering estimates and tabular data from 
standard handbooks. 

The primary items of interest are heat transfer coeffi- 
cients and initial conditions such as steam flow rates, inlet 
and outlet temperatures, circulation rates, etc. While 
some of these were directly obtainable from the boiler 
Technical Manual, others such as tube wall temperatures 
and riser exit quality were not. Therefore, a program was 
developed to permit the user to enter known data elements 
from the Technical Manual and other sources, and to obtain 
as an output the necessary values for the constants and 
initial conditions to ensure static balance at the speci- 
fied operating point. | 

The basic premise under which the program operates 
is that the model equations with all derivatives set to 
zero form a determinant set with the addition of several 
relations which apply under equilibrium conditions. These 
relations are: 

a. Riser and downcomer mass flow rates in : 


each circuit are equal. 


oy 





b. Total steam flow from the boiler equals 
the incoming feed flow rate. 

c. Total heat transfer to the incoming feed 
equals the heat transfer to the screen and the 
main bank, and the heat removed from the steam 
passing through the desuperheater. 

dad. The net mass transfer at the steam liquid 


interface in the drum is zero. 


B. PROGRAM EQUATIONS 


The total heat input to the boiler was given by 


Ee 2 Heat Release (Bw /sec- ft?) » Furnace Volowse (443) 


If the sensible heats of the oil and air are lumped with 


the heating value of the fuel, the heating value is 


Ho = Es /Ge 


The heat abSorbed by the furnace screen at steady state 


was calculated from 


Eai= Furnace Hear Apsonstin (Btu / sec - 4?) 
x Radiant Heat Absorbing Ourface (ft?) . 
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Since the flue gas temperature and the screen tube wall 
temperatures were provided, the heat transfer coefficient 
for radiant heat transfer to the screen could be determined 


from 


Yee Ena / { ( ( Tes +40 +0 ) = e (ae)"} 


~ 400 


This result holds for the conditions reported in the 
manufacturers technical manual. Degredation due to age or 
fouling must be accounted for through the use of appropriate 
correlations. 

With the enthalpy of the feed water at the economizer 
outlet and the enthalpy of the steam leaving the dry pipe, 
and the change in enthalpy of the steam flowing through the 
desuperheater, the heat transfer to the main bank was 


calculated from 


Ees = Coa (Hor Hwy ) — Era — Gor CHeyeur — Hos our) ; 


The metal temperature of the main bank riSers were obtained 


by calculating the constant 


kway ___Ame 


L182 % 10° 


and solving the heat transfer equation 
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We 2 13 
“Tm = { Ea3 / Kwy Poa : t +104 


The specific heat of the flue gas passing the main bank 


was 


Ca, = Eas / Ge Tr; - Tey) 


and the unknown constant in the gas side heat transfer 


relation was determined from 


Kes = Ees / Gee (Tan3 -Tra) 


In a Similar manner, the specific heat of the furnace 


flue gas is 


Cer = CEg-fat) / Gari 
and the water side heat transfer coefficient is 


a 4/3 3 
Kuz = Eas / Poa! (Tm4 —Toz ) 
The heat transferred in the desuperheater was: 


Eo2 = Gor CWster — Hpsour ) 
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With this value and the temperature difference the specific 


heat of the steam was computed from 
Cor >= Eon / Gor (Tor -T3 ) 


The heat transfer coefficient was calculated from the 


Dittus-Boelter correlation [ll] 


-o 
Kor= — = (.023) ea 5 oe, p.* 


The log mean temperature difference was 


»8 
LMTD = For /icor Gor 


and the metal temperature was determined by 


“Tas = (Tor 7 QToa ) / di-@) 
where: 


Q = r=Xp { Lmvd /(Tr-T03) 
With the metal temperature, the water side heat transfer 


constant was obtained from 


Kwr2 = Enz / CTms — Tw ) 
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Analysis of the Economizer is more complicated since it is 
a multipass heat exchanger with finned tubes. Additionally, 
since we are concerned with energy storage within the metal 
of the economizer, the internal and external heat transfer 
must be considered separately. 

The heat transferred was given by the change in enthalpy 
of the feed water, l.e., 


Ewt = Goa (Hwy — Hwa ) 


The specific heat of the water was 
Cw = Bwt / Goa (Tisy-Tes) 


As in the case of the desuperheater, the inside heat 
transfer coefficient was calculated from the Dittus- 


Boelter correlation as 


= 
Kwa = = (.023) (_t) Fr ° A 


From this the LMTD was obtained as 


LMTDd = Eur / Kw Gou 
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Finally the metal temperature was determined from 


Vg = (Ted - @tiws ) / Ci-A) 


B= €xp { LMT I'D / (Two Ear at 


The outside heat transfer constant was 


Kay = Ewt / Ge” (Tras -Twa) 


where 


The specific heat of the flue gas passing the economizer 


was 


Ces = Ews. / Ge CTry -Trs) 


The superheater heat transfer was obtained from the change 


in enthalpy of the steam, l.e., 


Foa= Got C Usn., - Hos ) 


With this value and the known metal and gas temperatures, 


the specific heat of steam, the specific heat of the flue 
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gas passing the superheater, and the inside and outside 


heat transfer constants were determined as 


Coa = Est / Gor CTor - Tos) 


Ceta = En / Ga (Tez - Tra) 


toa = Bod / Goa'® UmTy 
R22 = Epa / Gr’ (Tea - Ty) 
where: 
Tor - “Tha 
LMTO = “Tot 
| 
Tey - Tox 
and 


Teast = CTar + Tr) /iz.0 


The superheater and desuperheater pressure drop constants 


were evaluated from the known equilibrium pressures as 


Cy 
Cq 


C Poa -Por ) ose /Go4 
CPr2- Pos) Praz /6o2° 
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Similarly, the constants used in the calculation of the 


average steam densities were obtained from: 


Ca = CPoa+Por) /( A, * Tor ) 


and 


Cla 
\ = 
© CP», + Piy3 \ / (Pir +Bs) 
Finally the steam value flow constant was computed from 


C3 = Goy / (Poz Yewr ) 


C. CIRCULATION BALANCE 

Since they are not normally provided by the manufacturer, 
values for the friction factors, and bend and entrance 
loss coefficients must be assumed. The values do not appear 
to be particularly critical, and tables found in engineering 
handbooks may be used to determine reasonable values for 
these parameters. 

An iterative procedure must be used to adjust the various 
flow values involved to achieve a Steady state balance. 
Since the riser exit quality for boilers of this type 
generally lies in the range of .02 to .04, a Starting trial 


value for the circulation flow rate may be obtained from: 
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Ey = Xe Gy H 54 


where Xe is the exit quality, EL is the heat transfer to 


the riser, Gy is the flow rate, and HFG is the heat of 


Vaporization. With this trial value, the following 


balancing procedure is used: 
1. Set HW2 (Downcomer Enthalpy) = HW1l 
2. Calculate HWG = HW2 + ED2/G¢ 


3. Calculate the riser exit quality: 


Xo = Eas / CUgq + Hwa He ) Go 


4. Determine the non-boiling height: 


Waea = Axe CHwi-Hwe) /CHwt + Xo Hfg — Hw.) 


5. Compute the average riser density: 


1 § UxarHuse Ko Vég 
ran e { Xo Veq In ( VE t 1) a fs Hwee { 


6. Obtain the outlet density: 


Ke = 1.0 / CVE 4+ Xo Veq ) 
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10. 


a ae 


ae 


Compute the friction factor: 


eee ” 
* 1d = 2 log Cee | + Wworg.ss (32> 


Next the riser loss factor is computed: 


Ya = (Fae i) se aaeRe +Z3/aArkw + Zy J2Ar Ro 


Finally we may calculate Gg': 


Go = if ( Ru — Ras ) Hxe 2) / CCioo + v1) = 


where C100 is the downcomer loss factor. 


Compare Gg and Gg' if equal go on, otherwise 


update G@: 


and return to step 2. 
Assume a trial value for Gs (first pass only). 


Compute Xs: 


Xs = (Goi — Go Xe) /&s 
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13-19 Same as for the main bank risers. Return to 


12 if reiteration necessary. 


20. Obtain HW2' from: 


/ Gw2 
Wir = Goa Ge CHa Hus) + wa 
21. If HW2' = HW2 go on, otherwise update HW2: 


4 
Hw. —Hw2 


+ Hw2 
Z 


Hw2r = 


and return to step 2. 


Upon exit from the iteration loop, equilibrium values 

of the flow parameters will have been established. With 
these and the previously calculated heat transfer parameters, 
the program generates the list of constants and integration 
initial conditions required by the dynamic model. The 


program listing and a sample output 1S given in Appendix B. 





V. RESULTS AND CONCLUSIONS 


A. OPEN LOOP RESPONSES 

The output curves for 10% increases and decreases in 
throttle valve opening and fuel flow are given in Figures 6- 
33. Ramp inputs of approximately 1.5% per second were used 
to simulate flexibility test conditions. Only one input 
parameter was varied for each run with all others held 
constant. The reSultsS were similar to those presented in 
[3] and [4] insofar as general trends were concerned. Since 
the boilers involved were of widely different geometries 
and sizes, a quantitative comparison was not possible. 
Since CSMP uses only Single precision arithmetic packages 
and integration algorithms, the high frequency noise in the 
outputs 1S attributable to the roundoff error inherent in 
the wide variations in magnitudes of numbers used in the 
model, e.g., quality (order of 0.01) and riser mass flow 
rate (order of 1000). In common with the results of [3] 
and [4], the "swell effect" in response to an increase in 
throttle opening is negligible. It should be noted that the 
swell exhibited in the data of Whalley [3] 1S more apparent 
than real, inasmuch as the scales in his graphs are very 
magnified and the rise is actually only 0.01 inches. The 
Sudden increase in riser mass flow rate shown in [3] is 


also not present in the current model. 
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Since "Shrink" and "swell" effects are of some concern 
to controls designers, an attempt was made to model for 
Such effects. Four possible causes for the phenomenon 


were considered. 


1. Homogeneous Nucleation in the Steam Drum Liquid 


This effect is caused by a difference between the 
Saturation temperature corresponding to the actual drum 
pressure, and the actual liquid temperature. The super- 
heat required for such nucleation iS quite high however and 
nucleation will occur at the heating surfaces at much lower 
superheats. Therefore this was considered to be an unlikely 
cause for the effect. 

2. NucleateBoiling at the Drum Wall-liquid Interface 

The driving force for this phenomenon is also the 
temperature difference between the drum surface and the 
Saturation temperature. Since the temperature difference 
will be slight except under large drum pressure excursions, 
the impact of this factor should be slight. The Rohsenow 
correlation [5] may be solved in terms of a rate of vapor 


generation: to yield 


Cp, AT : - Tv) 


Bheg = 1 Ga tg Pep Ae o 


It was assumed that the vapor formed would be uniformly 


distributed bubbles which would rise to the surface with a 
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velocity given by [14] 


U, = 3V 29 Dd 


which is a conservative estimate since interaction effects 
tend to increase the rise velocity. The bubble diameter 


is given by [14] 


D= -24, P'S 


This gives UL as .25 ft/sec. The average liquid depth is 

15 inches, so the transit time is =~ 5 seconds. The effect 
was modeled crudely by assuming the entrained bubble 

mass was 5.0(q/hfg) lbm. The impact of this model on shrink 
and swell values was negligible, since no large pressure 
excursions were developed. 


3. Expansion and Contraction of Varpor Bubbles Entrained 
in the Riser Fluid 


An attempt was made to model this effect by including 
the riser tubes in the control volume used as the basis for 


liquid level determination. The equation used was 





Vo. = Mwi /0 w + (fe -1)\re + (Kae —1) be 


where Ralx and Ra2x are the average riser densities computed 


using the specific volume for the saturated vapor corresponding 
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to drum pressure, while Ral and Ra2 were computed from 
the liquid temperature saturation conditions. The 
effects were inSignificant. 


4. Transient Deviation of Riser and Downcomer Mass 
Flow Rates 


The small swell effect obtained by Whalley [3] 
can be attributed to the temporary increase of riser mass 
flow rate over downcomer mass flow rate. The period over 
which this imbalance could be sustained is short and would 
fail to account for the longer duration of "shrink" and 
"swell" effects observed in practice. It may, however 
trigger one of the two phase flow instabilities common to 
natural circulation loops e.g., flow oscillations [10,13]. 
The analysis of these effects was not exhaustive, and the 
modeling attempts were crude. In the case of the first two, 
however, the contribution to swell is probably not signifi- 
cant. Modeling difficulties with the latter two effects 


precluded any specific conclusions. 


B. CONCLUSIONS 

As with all simulations, the information available from 
the model is extensive. Notably, in this study, main bank 
and water~-Screen circulation values can be separately 
evaluated. Hence, the contribution of each of these ele- 
ments to overall boiler operational response and efficiency 
can be determined. 

The model simulates the responses of a Naval boiler 


quite well. Unfortunately, open loop data is unavailable 
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for the actual boiler and insufficient information was 
available to develop a model of the existing LHA-1 control 
system as a basis for comparison studies. As such, however, 
the model forms a basis upon which future studies involving 
modern control system designs can be established. Moreover, 
while considerable refinement of the model is possible, 

the advantages of greater accuracy must always be weighed 
against increased computation costs. In this regard, the 
implementation of more efficient integration algorithms 


would be particularly beneficial. 
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